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Yang-Mills theories in 2+1 (or 3) dimensions are interesting as nontrivial gauge theories in their 
own right and as effective theories of QCD at high temperatures. I shall review the basics of our 
Hamiltonian approach to this theory, emphasizing symmetries with a short update on its status. 
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We will show that the calculation of the vacuum wave function for Yang-Mills theory in 2+1 
qj dimensions is in the lowest order of a systematic expansion. Expectation values of observables 

7^ , can be calculated using an effective interacting chiral boson theory, which also leads to a natural 

expansion as a double series in the coupling constant (to be interpreted within a resummed per- 
turbation series) and a particular kinematical factor. The calculation of the first set of corrections 
in this expansion shows that the string tension is modified by about —0.3% to —2.8% compared 

J^-j to the lowest order value. This is in good agreement with lattice estimates. 
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1. Introduction 

As many of you know, Yang-Mills theory in 2 + 1 dimensions has been the focus of research 
as well as an object of fascination for my collaborators and myself for many years. In the first part 
of this talk, I shall give an update on the status of our attempts to understand this theory. In the 
second part, I shall zero in on our recent calculations of corrections to the string tension. This latter 
part is more than just a calculation of specific results; it also formulates a systematic calculational 
framework which may prove to be a gateway to newer results such as glueball masses. 

Let me begin by recalling why this theory is of interest to us, although for a QCD-savvy audi- 
ence such as this, such a reminder is perhaps redundant. Yang-Mills theories have, of course, been 
a great puzzle for many decades. Since the realistic case of 3 + 1 dimensions is highly nontrivial 
and difficult, physicists have often turned to the time-honored path of seeking guidance from lower 
dimensions. Yang-Mills theories in 1 + 1 dimensions are exactly solvable, but they have no propa- 
gating degrees of freedom and, as a consequence, are a little too trivial for us. In 2 + 1 dimensions, 
there are propagating degrees of freedom, the theory has nontrivial dynamics, yet seems to be 
within reach of some level of analytical investigation. This is aided by the fact that the theory has 
a dimensional coupling constant and, as a related fact, because it is super-renormalizable. There 
is also an interesting physical context in which Yang-Mills theory in 2 + 1 dimensions, in its own 
right, is useful. One can use Yang-Mills theory in 3 Euclidean dimensions as the high temperature 
approximation to real QCD; calculations in the Minkowskian 2+1 dimensional theory can then be 
continued to the Euclidean signature to make predictions about the magnetic screening mass and 
other related phenomena. 

2. A short review and status report 

Our approach to this problem has been to use Hamiltonian analysis and to solve the Schrodinger 
equation for the vacuum wave function [1,2]. Generally, Hamiltonian methods are not easy to im- 
plement in a field theory with its infinite number of degrees of freedom and regularization and 
renormalization issues. But for this particular problem, it has some advantages, because, at con- 
stant time, we are dealing with two-dimensional gauge fields and, for these, our experience with 
conformal field theory is a useful guide. 

As is standard in a Hamiltonian analysis, we then start with the Ao = gauge. The spatial 
components of the gauge potential can be combined as A = A z = \ (A \ + /A 2 ) , A = A| = \ (A \ — /A 2 ) . 
The crucial first step for us is the parametrization of these potentials as 

A = -dMM~\ A=M^ l dM\ (2.1) 

Here M is an SL(N, C) -matrix for an SU (N) -gauge theory. More generally, M is a complex matrix 
which is an element of G c , the complexification of G which is the Lie group in which the gauge 
transformations take values. 

Time-independent gauge transformations act on M via , 

M{x) — ► g{x)M{x) , g(x) E SU(N) (2.2) 
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The hermitian matrix H = M^M is thus gauge-invariant. The fact that the gauge transformations 
are made homogeneous, as in (2.2) is why the parametrization (2.1) is so useful. 

Wave functions ^(H) are gauge-invariant and depend only on H = Af^M; this may be thought 
of as the statement of the Gauss law. Further, the Jacobian of the transformation A, A ^ H can be 
explicitly evaluated and the inner product for the wave functions can be obtained as 

(1|2) = J dn(H) exp(2c A S wzw (H)) *F*¥ 2 (2.3) 

where S WZW (H) is the Wess-Zumino-Witten action for the hermitian field H and dji(H) is the Haar 
measure for H viewed as an element of SL(N, C)/SU (N). The number ca is the quadratic Casimir 
value for the adjoint representation; it is defined in terms of the structure constants f abc (of the Lie 
algebra) as c A 8 ab = f amn f bmn and is equal to N for SU{N). The WZW action is given by 

= ^J TrtfHdH-^ + ^J e^ a Tv{H- l d^H H~ l d v H H~ l d a H) (2.4) 

Actually we can strengthen the statements given above a little further and argue that the wave 
functions *P and the Hamiltonian can be taken to be functions of a scaled version of the current of 
the WZW action, namely, of 

2 

J=-dHH~ l (2.5) 
e 

One way to see this is to consider the Wilson loop variable, 

W(C) = Tr &>e~k A = Tr ^exp (^f J> ) ( 2 - 6 ) 

where, in the second step, we have simply used the parametrization (2.1). Since all gauge-invariant 
quantities can be constructed from the Wilson loop variable (by suitable choices of the contour of 
integration C), we may restrict attention to functions of /. 

Taking the wave functions to be functions of the current /, functional derivatives with respect 
to A, A may be obtained via the chain rule of differentiation. This leads to the expression for the 
Hamiltonian operator in terms of the / as = + J(f[ with 

j% = m [j a ®*^ + l[ 1 , 2 - * * +^ [:dJ a (z)dJ a (z): (2.7) 
J oJ a (z) 7i J w , z (z-w) z oJ a (w) oJ a (z) 2 J z 

J c {w) 8 8 



= ie ( fabc—? 



r(z- w) 8J a (w) 8J b {z) 

where m = e 2 CA/2n = e 2 N/2n. Notice that the first term of the Hamiltonian is counting powers of 
/ in the wave function on which it acts, giving an energy contribution m for each /. This will turn 
out to be the essence of the mass gap in the theory. 

An aside on regularization 

The Hamiltonian operator involves products of /'s, 8/8J's and nonlocal functions such as 
(z — w) _1 . As in any field theory, all calculations have to be done with proper regularization. The 
one we have used involves a thickening of the Dirac 8 -function, 

S^Xu.w) cr(w,w;e) = ^exp ( J u ~^ \ (2.8) 
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This leads to the replacement of the unregulated Green's function for d by its regularized version 
%ia(x,y)\ i.e., 

G(x,y) = — r - — r — > ^ma(x,y) [ G(x,u)a(u,y;e) H(u,y)H~ l (y,y) 

* i^[^- e ~^* le ( H ( x &r x m)^ (2.9) 

The parameter controlling the regularization, £, acts as a short-distance cut-off. Expression (2.9) is 
to be used on functionals where the point- separation of various factors is much larger than y/e. 

While the parametrization (2.1) has many advantages, it has an ambiguity in that (M, M^) 
and (MV(z), V(z)M^) give the same gauge potentials A, A, where V, V are, respectively, anti- 
holomorphic and holomorphic in the complex coordinates z = x\ + 1x2 and z = x\ — 1x2. To avoid 
this ambiguity of parametrization, physical observables in the theory should have invariance under 
H — > V(z)HV(z). The regularization used in (2.9) respects this "holomorphic invariance"; in fact, 
factors like (H(x,y)H~ l (y,y)) appear in (2.9) precisely for this reason. 

All the results we present here have been checked using regularized expressions, with a single 
regulator from beginning to end, to avoid possible conflicts between regularizations. The story of 
regularization, at the level of calculations, is more involved than what we have indicated here; for 
more details, see [1, 2, 3]. 

One interesting result can be read off from what we have done so far. By comparison with 
resummed perturbation theory, we can see that m serves as the propagator mass for gluons (which 
are generated by J). This can then be identified with the magnetic screening mass if we interpret the 
theory as the high temperature limit of QCD. For the gauge group SU (2), m = e 2 / 71 « 0.32<? 2 . There 
are a number of other calculations for the magnetic screening mass, so it is interesting to compare 
these values. Table 1 shows that while there is some variation, the values are not outrageously 
different from e 2 / 71. 

In our earlier work [1, 2], we also solved the Schrodinger equation Jif^o = for the vacuum 
wave function ^0 to the leading order in a strong-coupling expansion to obtain ^0 = 2 » where 



Table 1 


. Comparison of magnetic mass calculations 


m/e 2 


Method 


0.25 


Resummation of perturbation theory [4] 


0.35 


Lattice, common factor for glueball masses [5] 


0.51 


Lattice, maximal abelian gauge [6] 


0.52 


Lattice, Landau gauge [6] 


0.44 


Lattice, A3 = 2 gauge [6] 


0.38 


Resummation of perturbation theory [7] 


0.28 


Resummation of perturbation theory [8, 9] 


0.37 


Gauge-invariant lattice definition [10] 


0.32 


Calculation via our Hamiltonian method 



4 



Yang-Mills (2+1) 



V. P. NAIR 



S(H) = J dJ a 



1 



(m + \fn 



/m--V 2 ) 

Notice that if one restricts to modes of / with momentum <C e 2 , 



dJ a + 0(J 3 



S(H) 



[ dJ a dJ a = / d 2 x FfiFfi 
2m J 4g 2 J lJ lJ 



(2.10) 



(2.11) 



Thus the computation of expectation values reduces in this limit to a calculation in a Euclidean two- 
dimensional Yang-Mills theory with a coupling g 2 = me 2 . The expectation value of the Wilson loop 
variable for the representation R then takes the form 



(W(C))~exp(-a^c) 
where srfc is the area enclosed by the curve C. The string tension Or is given by 



caCr 
An 



(2.12) 



(2.13) 



This compares favorably, to within 1 — 3%, with the latest lattice calculations [11, 12], see table 2. 

I have not described in detail the process of solving the Schrodinger equation; this is be- 
cause there is another argument which can be used to obtain the results to this order [2]. For 
this purpose, absorb the factor oxp(2caS wzw (H )) from the inner product into the wave function by 
\j/ _ e -c A s wzw {H)^^ Then the Hamiltonian acting on <I> is given by 



Table 2. Comparison of predictions from (2.13) with lattice calcuations 



(2.14) 



Group 


Representations 




k=l 


k=2 


k=3 


k=2 


k=3 


k=3 




Fund. 


antisym 


antisym 


sym 


sym 


mixed 


SU(2) 


0.345 
0.335 












SU(3) 


0.564 
0.553 












SU(4) 


0.772 
0.759 


0.891 
0.883 




1.196 
1.110 






SU(5) 


0.977 
0.966 












SU(6) 


1.180 


1.493 


1.583 


1.784 


2.318 


1.985 




1.167 


1.484 


1.569 


1.727 


2.251 


1.921 


SU{N) 


0.1995 N 














0.1976 













The upper entries follow from (2.13), the lower entries are the lattice values, k is the rank of the 

representation. 
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We shall now consider this expression when H can be approximated as H = e ta(pa « 1 + t a (p a + • 
this small <p limit is appropriate for perturbation theory. The new Hamiltonian is then 



if 



8(j> a 8<l) a 



+ f (m 2 -V 2 )(j) a 



(2.15) 



where </> fl = y/c A (-V 2 )/8nm (p a . We see that the leading term in corresponds to a free field 
of mass m (actually dim G fields, counting the multiplicity due to the index a.) The field may 
be taken as representing the gluon (actually a gauge-invariant version with some resummations 
involved) and this is consistent with what we said about the gluon mass being m. 

We may now trivially calculate the vacuum wave function in this approximation as 



<&o ~ exp 



If this is transformed back to the ^-version, we find 



« exp 



Km J 



m- 



1 



-V 2 <\) a 



Wm 2 -V 2 



m 



{dd<p a )+-- 



(2.16) 



(2.17) 



This was obtained in an approximation of small <p. However, we do know that the full wave 
function must be a function of the current /. So we can ask the question: What function of / 
can we write such that it reduces to (2.17) in the small (p approximation when / « (2/e)d<p7 
The answer, to quadratic order in /, is clearly the solution (2.10). Notice that this argument only 
depends on the form of integration measure which, in turn, is determined by a two-dimenional 
anomaly calculation. So it is rather robust. Of course, a regularized theory is required so that we 
can carry out the transformations *P <-> <& without worrying about additional effects coming up. 

Another interesting line of development has been about glueball masses. The authors of [13] 
have followed the general line of reasoning we have outlined, but used a somewhat different wave 
function given by 



exp 



-i 

Am J 



dJ K[L] dJ 



K[L} = 



Ji{aVl) 

VLJi(4y/L) 



(2.18) 



where L = $>d jm 1 and Ji, J2 are Bessel functions of orders 1 and 2 respectively. The kernel K is, 
despite appearances, very close to the kernel in (2.10), as is clear from the graphical comparison 
in figure 1 (I thank A. Yelnikov for this picture). Glueball masses were obtained by calculating 
the two-point function for different color-singlet composite operators, characterized by spin, parity 
and charge conjugation properties (j PC -notation). The results, in units of ^ S /Gf, are shown in table 
3. Again, there is reasonable agreement with the lattice data of reference [11, 14]. 

There have been a number of other developments as well. The formalism has been extended to 
the Yang-Mills-Chern-Simons theory and issues related to dynamical mass generation and screen- 
ing of Wilson loops have been analyzed [15]. 

Another interesting question has been about the screening of Z/v-invariant representations. The 
energy of a glue-lump state made of a heavy scalar field and the gluons can be estimated in the limit 
of large e 2 and large N [3]. This gives the energy at which an adjoint string breaks as E* « 7.92 m, 
for SU (2), to be compared with the lattice value of E* « 8.68 m [16]; we are off by about 8.8%. 
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Figure 1. Comparison of KKN and LMY kernels for the Gaussian term of the vacuum wave 
function as a function of the momentum p 

Considering the difficulties in calculating both in our method and the lattice, this may be taken 
as not too bad. The calculation is more a qualitative demonstration of feasibility: it is possible to 
incorporate screening in our formalism. 

We have also worked out the extension of the formalism to the case of space being S 2 , i.e., 
for R x S 2 [17], as a prelude to calculations on the torus which may be used to analyze finite- 
temperature effects. In this regard, I also mention the recent paper by Abe, where some of the torus 
formalism is developed and arguments on deconfinement are presented [18]. 



3. Corrections to string tension 

In the second half of the talk, I shall focus on corrections to the string tension formula (2.13) 

Table 3. Comparison of glueball mass estimates from (2.18) and lattice calculations. 



State 


LMY Calculation 


Lattice 


0++ 


4.098 


4.065 ±0.055 


o++* 


5.407 


6.18±0.13 


0++** 


6.716 


7.99 ±0.22 




7.994 


9.44 ±0.38 


0~ 


6.15 


5.91 ±0.25 


Q— * 


7.46 


7.63 ±0.37 


Q ** 


8.77 


8.96 ±0.65 


2++ 


6.72 


6.88±0.16 


2 ++* 


7.99 


8.62 ±0.38 


2^-\ — h** 


9.26 


9.22 ±0.32 


2+" 


8.76 


8.04 ±0.50 


2~ 


8.76 


7.89 ±0.35 


2 +-* 


10.04 


9.97 ±0.91 


2 — * 


10.04 


9.46 ±0.46 
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[19]. This is important for several reasons. First, we want to see the vacuum wave function (2.10) 
and the corresponding string tension (2.13) as the lowest order in a systematic expansion scheme. 
Secondly, notice that, in the large N limit, our formula differs from the lattice value by about 1%. A 
different lattice calculation shows the difference to be about 1.55% [20]. In both cases, the lattice 
calculations are considered accurate enough that the deviations are statistically significant. So we 
need a systematic method which has the potential to calculate corrections. 

There are two types of corrections possible as exemplified by the diagrams in figure 2. For the 
first set, we can calculate the corrections to the propagators (or the coefficient of dJ dJ in the wave 
function) and then use the result in the calculation of the Wilson loop. This set of corrections will 
be independent of the representation R of the Wilson line. The second set would involve vertex 
type corrections and can be sensitive to the representation of the Wilson line. We will leave out 
these for now, focusing on the first set. (But note that the second type are related to screening and 
string breaking problems, which we have commented on already. Also, since these /^-dependent 
corrections are left out, comparison with lattice analysis is best made at large N.) 

The corrections are expected to be purely numerical, so, at first glance, seem to be impossible 
to control in any way. However, we shall see that a systematic expansion is possible along the 
following lines. 
Step 1 

We rewrite the derivation of the vacuum wave function as a recursive procedure for the solution 
of the Schrodinger equation from which it will be clear that (2.10) is the lowest order result in a 
systematic expansion. For this purpose, we will treat m and e as independent parameters, setting 
m = e 2 CA/2n only at the end of all calculations. For ^q^Fo = e F , this gives F as a power series in 
e, of the form 

2 

+ J Aa 3 a 4 (*l,X2,X3,X 4 ) (x l )j a '(x 2 )j a ^X 3 )j a '(x 4 ) + ■■■ (3.1) 

p 2 \ etc., are determined recursively in powers of e, with each function having an expansion 

fa}a 2 {xi,X 2 ) = fol a2 (x U X 2 ) + <? 2 f[ \ a2 (x 1 ,X 2 ) +••• 

fola 2 a 3 (x\,x 2 ,x 3 ) = £ } ai a 2 a 3 (*i,x 2 ,x 3 ) + fi!^ (x\,x 2 ,x 3 ) +■■-, etc. (3.2) 
To the lowest order, with q,q denoting the momenta, 

/(fiajC*!^) = <5a,a 2 [-f /( m + E q)] Xl ,* 2 ~ S a,a 2 f /2m (3.3) 




Figure 2. Corrections to Wilson loop (rectangle); wavy lines are (//) propagators. 
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(i\ faia 2 a 3 ift 



+ cycl. perm. 



(3.4) 



which leads to (2.10). To calculate corrections to order e 1 , we need the lowest order results for 
and /( 4 ) . For example, 

hh(h ~h) 
(m + E kl )(m + E kl ) 

There is a similar, but more involved, expression for /W, see [19]. 
Step 2 

We are interested in corrections to since that determines the string tension. From the recursive 
solution of the Schrodinger equation, we find 

j2l /1 h \ ?2 



e 1 



(3.5) 

Seemingly, this is a 113% correction, but, as we shall see, other "loop" corrections are important. 

We want to emphasize that, even though we use recursion in e, our procedure is different from 
perturbation theory since m is included exactly in the lowest order result for F\ more on this later. 
Step 3 

Since the measure of integration has the WZW action, we can transform the functional integration 
over ^o^o = e F into the integration over a two-dimensional chiral boson field <p, <p (not the <p we 
used in parametrizing H as e tCl(pa ). In other words, for an observable &, 



= J dn(H) e 2cASw ^ H) e F{J) 0{J) = J [dcpdcp] e~ s ^ 0(^2n/mc A fpt a (p) (3.6) 

where the action S((p) is given by 

S((p) = J (Z 2 (pd(p + Zi(f)C(p) - F{Ziyfl7i/mc A q>t a (p) (3.7) 

(Zi ,Z2 are renormalization constants, we will not discuss these further here; see [19].) 1 Notice that, 
effectively, the current J a is replaced by Z\ ^2%/mcA q>t a (p. The function F{Z\ ^2n/mcA q*t a (p) 
contains vertices, F^ with two currents (quartic in <p, <p), F^ with three currents, etc. For exam- 
ple, we may diagrammatically represent F^ 2 \ with two <p's and two (p's, as 

Fi2) = ^J(<pt a <p) x fV{x,y)(<pt«<p) y = 



The corrections to F^ 2 \ which is what we are interested in, may be viewed as Wilsonian renormal- 
ization corrections to the quartic vertex in this two-dimensional field theory. 
Step 4 

In computing the corrections to F^ 2 \ we can treat the vertices F^ 3 \ F^\ etc., perturbatively since 
they carry powers of e. However, the lowest term in the vertex F^ 2 \ corresponding to f^ iai (x\ , JC2), 
has no powers of e and hence its contributions must be summed up. This means that the current- 
current correlator becomes 

w*w |r'*"-i(t) (38) 

l C may be set to zero at this point, it was included to show how Z\, Z2 appear in S(<p). Also we may note that the 
transformation to <p, (p is analogous to the fermionization of the WZW model. 
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{<pt a tp(x) <pt b <p(y)) = x<^ J>y 




Figure 3. The current-current correlator including all contributions from F^ 






Figure 4. Corrections from F^ summed up as a factor of m/Ek (shaded circle at vertex) and 

sample renormalization diagrams 



Here Ek = Vk 2 + m 2 ; the m/Ek factor arises from the summation of corrections due to Fq 2 \ shown 
diagrammatically in figure 3. Our computational strategy may then be summarized as follows: 

1. Construct loop diagrams generated by F^ (3 factors of (pt a q>) and F^ (4 factors of (pt a q>). 

2. They can have arbitrary insertions of F^'s, leading to a factor of m/Ek, as in figure 4. 

3. There are also renormalizations (due to F^) we have to take into account. 

4. Sum up F^ insertions in all diagrams (of order e 2 ) generated by F^ and F^ 4 ). 

5. Classify and group these by the number of factors of m/Ek. 

We will compute corrections to order e 2 and up to 4 powers of m/Ek. Denoting the factors 
of m/Ek by shaded circles at the vertices, the corrections to the low momentum limit of may 
summarized as in figure 5. We show the coefficients of q 2 /2m, for small q, q. 

Let C n denote the partial sum of corrections up to terms with (m/Ek) n , starting with Co = 
1.1308 from the recursive procedure (3.5). We then find Ci = 0.5496, C 2 = 0.2730, C 3 = 0.0373. 
There is a small ambiguity in an integral for the last diagram in figure 5 [19], so that C4 is in the 
range —0.05843 to —0.00583. The partial sums are systematically decreasing in value, showing 
that the ordering of diagrams by powers of m/Ek does constitute a sensible expansion. Also, the 
corrections to the order we have calculated are small. For the string tension, we find 



caCr 
4n 



0.02799 + ■ 
-0.0029 + •• 



) 



(3.9) 



This correction, of the order of —2.8% to —0.03%, is entirely consistent with lattice calculations. 
Terms of order (m/Ek) 5 are expected to contribute at the level of a fraction of 1%, and likewise for 
diagrams with two or more current loops. 



10 



Yang-Mills (2+1) 



V. P. NAIR 




-0.58118 -0.47835 0.20169 -0.23569 




0.02083 -0.06893 




-0.01216 0.06824 (-0.1037) to (-0.166) 

Figure 5. Corrections to the low momentum limit of the vertex 



4. Discussion 

Perhaps the most important point to be highlighted is the fact that one has a systematic expan- 
sion and a scheme for calculations. At first glance, the possibility of any expansion scheme in 2+1 
dimensional Yang-Mills theory seems very remote, to say the least. The coupling constant e 2 does 
not constitute an expansion parameter. It simply tells us that modes of momenta <C e 2 should be 
treated nonperturbatively, while modes of momenta » e 2 can be treated perturbatively; it is at best 
a marker for this separation. Nevertheless, from what we have discussed, a systematic calculational 
procedure is possible. This is done by keeping m and e 2 as independent until the end, giving a para- 
metric way to classify various contributions to any calculation. Further, there is a kinematic factor, 
of the form (m/Ek) n , which helps in the further classification of corrections. The corrections to the 
lowest order results for the string tension then come out to be rather small. 

Resummation calculations have also indicated some sort of smallness for the corrections. The 
magnetic mass results quoted came from calculations to "one-loop order" for the Schwinger-Dyson 
equations. In other words, a particular infinite set of diagrams were summed up which corre- 
sponded to one-loop order viewed as a skeletal diagram. A "two-loop" extension of this, with no 
obvious reason, gave a reasonably small correction [21]. Notice also that the values from resum- 
mation (coming from the perturbative side) and our value (coming from the large e 2 side) are not 
very different. There is some sort of stability, almost "unreasonable stability", to the calculations in 
Yang-Mills (2+1). This suggests that there is a more natural scheme of calculation for this theory, 
some features of which we are discovering, perhaps the hard way. 

Thinking in the language of ordinary perturbation theory, we may note that our analysis in- 
volves at least three different resummations. First, the definition of H in terms of A, A is an infinite 
series with nonlocal terms. Second, in including m in there is another resummation, since m 
is really of order e 2 . Thirdly, we also carry out the summation of contributions due to F^ to get 
factors of mjE^. An a priori attempt at these resummations, even with some guidance from the 
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Hamiltonian analysis, is clearly very difficult. 

Where do we go from here? Higher order corrections to the string tension, a better computation 
of glueball masses, extension to the torus and finite-temperature effects as well as the inclusion of 
matter fields are some of the questions we could now explore systematically. 

I thank the organizers for inviting me to this very interesting workshop and for accommodating 
my difficult schedule. 
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